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ABSTRACT 

Discovery of the cosmic reionization epoch would represent a significant milestone in cosmology. We 
present Keck spectroscopy of the quasar SDSS 1044-0125, at z = 5.73. The spectrum shows a dramatic 
increase in the optical depth at observed wavelengths A > 7550 A, corresponding to z a b s > 5.2. Only a 
few small, narrow transmission regions are present in the spectrum beyond that point, and out to the 
redshifts where the quasar signal begins. We interpret this result as a signature of the trailing edge 
of the cosmic reionization epoch, which we estimate to occur around (z) ~ 6 (as indeed confirmed by 
subsequent observations by Becker et al. ), and extending down to z ~ 5.2. This behavior is expected 
in the modern theoretical models of the reionization era, which predict a patchy and gradual onset of 
reionization. The remaining transmission windows we see may correspond to the individual reionization 
bubbles (Strbmgren spheres) embedded in a still largely neutral intergalactic medium, intersected by the 
line of sight to the quasar. Future spectroscopic observations of quasars at comparable or larger redshifts 
will provide a more detailed insight into the structure and extent of the reionization era. 

Subject headings: cosmology: observations - galaxies: formation — galaxies: quasars: individual 
(SDSS 1044-0125) 



1. INTRODUCTION 

There has been a great progress over the past several 
years in our understanding of galaxy evolution and forma- 
tion. Sa mples of normal ga laxies are now studied out to 
z ~ 4.5 ( Steidel et al. 1999), and several g alaxies are now 
known at z > 5 (sec Stern fe Spinrad 1999] for a r e view and 
referenc es). Quasars at z > 5 (Stern et al. 2000, Zheng e\ 



seen. 



In this Letter we present evidence which suggests that 
we are already probing the trailing end of the reioniza- 
tion era, at z ~ 5.5 ± 0.3 or so. The evidence is based 
on the high S/N, Keck s pectroscopy of t he quasar SDSS 
1044-0125 discovered by Fan et al. 2000|. Subsequent re 



cent observations of quasars at z > 6 by Becker et al. 2001 



provide an additional evidence that the reionization era in- 
deed occurs around z ~ 6. Taken together, the available 
data support a picture of an extended and patchy reion- 
ization era, ending at z ~ 5 — 6. 

2. OBSERVATIONS AND DATA REDUCTIONS 

Our low-resolution spectra were obtained on the 
W. M. Keck Observatory 10-m telescope (Keck I) on 
UT 2000 December 30, u sing the Low Re solution Imag- 
ing Spectrometer (LRIS; Okc et al. 1995). The obser- 
vations were obtained with the 400 lines mm -1 grating 
(Ablaze = 8500 A) through 1.2 arcsec slitlets, with two 
different slitmasks, with a mean dispersion of w 1.86 A 
pixel -1 , and a GG495 long-pass order-sorting filter. The 
first set, totalling 2400 s of integration, was at a slit po- 
sition angle PA = 249.6°, covering the quasar spectrum 
~ 6240 A — 1 fim, and a mean airmass rj 1.08. The sec- 
ond, with 3600 s of integration, was at PA = 137° (very 
close to parallactic angle at the time), covering ~ 7300 
A - l /im, and a mean airmass 1.13. The differential 
slit losses are estimated to be negligible for our purpose. 

1 Based on the observations obtained at the W. M. Keck Observatory which is operated by the California Association for Research in 
Astronomy, a scientific partnership among California Institute of Technology, the University of California and the National Aeronautics and 
Space Administration. 
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al. 2000, Fan et al. 2000, and references therein) also repre- 
sent a valuable probe of both galaxy and structure forma- 
tion, and the intervening primordial intergalactic medium 
(IGM). 

The observational frontier is now shifting to the for- 
mation of the first objects, protogalaxies and primordial 
AGN, which is generally expected to occur some time in 
the redshift interval z ~ 6 — 15 or so. As the first sources of 
UV radiation turn on, they reionize the universe, ending 
the "dark ages" which start at the recombination epoch 
(z ~ 1100). In this "cosmic renaissance" (effectively, the 
start of the galaxy formation epoch) the universe under- 
goes a phase transition from being neutral to being mostly 
ionized. 

Detection of the reionization epoch would thus be a ma- 
jor cosmological milestone. The standard observational 
test is the prediction of an extended, optically thi ck ab- 
sorption due to n eutral hydrogen at X res t < 1216 A ( |Gunn 
fc Peterson 1965). A limit to this effect at z as 5 was pub- 



lished by 3ongaila et al. 1999. To date, only a gradual 
thickening of the absorption due to the Lya forest was 
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Reionization Epoch 



Individual integrations were dithered along the slit. Data 
were reduced in IRAF, using standard slit spectroscopy 
procedures. Ne+Ar arc lamp spectra obtained through 
the masks were used for wavelength calibrations, and 
the w avelength zcropoints adjusted using telluric emission 



mate a better redshift for the quasar: z — 5.73±0.01. This 



is less than 5.80 originally estimated by Fan et al. 2000 



for the followin g reason. The objec t was found to be a 
BAL quasar by [Vlaiolino et al. 200 1|. We believe that its 



lines. The night was photometric, but unfortunately no 
flux standards were observed with these slitmask+grating 
combinations, and we used an average of archival response 
curves for this grating obtained earlier. Our spectroscopic 
magn itudes are in an excellent agreement with the CCD 



Lya line is nearly completely absorbed, and that Fan ei 



al. 200C mistakenly interpreted the red half of the |N VJ 



photo xietry presented by Fan et al. 200C . 



The combined LRIS spectrum is shown in Fig. 1. It 
shows a dramatic drop due to the Lya absorption at 
A bs 8100 A, and a second discontinuity at A 6 S < 6900 
A, due to the Ly/3 forest (this may be the strongest detec- 
tion of the Ly/3 drop observed to date) . In order to esti- 
mate the possible continuum level in the absorbed region, 
we use three power-laws, f„ ~ v a , with a = [0, —0.5, —1], 
which span a plausible range for quasars. They have been 
normalized to /„ = 72 /iJy at 9000 A, and are shown as 
dashed lines in Fig. 1. 

Our high-resolution spectra were obtained on the 
W. M. Keck Observatory 10-m telescope (Keck II) on 
UT 2000 April 28, UT 2001 January 1 and 2, and UT 
2001 Marc h 24, using the E chelle Spectrograph and Im- 
ager (ESI; |Miller et al. 2000|) . We used the Echelle mode 
which contains ten orders, with a complete optical wave- 
length coverage, from ~ 3900 A to ~ 10900 A. The in- 
strument has a spectral resolution of 11.4 km s -1 pixel -1 , 
and a mean dispersion in the wavelength region of interest 
here of « 0.154 A pixel -1 . A total of 11 exposures of 1800 
s each were obtained, some in slightly non-photometric 
conditions. Data were reduced using standard procedures. 
We used the program MAKEE (written by T. Barlow) to 
reduce the spectra. Individual exposures from each night 
were combined prior to spectrum extraction using a rejec- 
tion algorithm to remove cosmic rays. Spectra were then 
optimally extracted. Exposures of bright stars were used 
to provide the spectrum traces (necessary due to a heavy 
absorption present in the quasar spectrum). Dispersion 
solutions were found from exposures of arc lamps, spectra 
were corrected to the Heliocentric system, and the wave- 
lengths were transformed to vacuum values. The spectra 
for each night were flux calibrated using a single response 
curve measured during one of the nights, and averaged 
using the exposure time weighting. 

Since the flux zero-points for the ESI data are uncertain, 
we convolved both ESI and LRIS spectra with Gaussians 
with a = 20 A, thus bringing them to effectively the same, 
very low resolution. From the ratio of these spectra we de- 
termined the flux correction factor, and applied it to the 
ESI data. 

The final ESI spectrum is shown in Fig. 2. The ab- 
sence of flux (save for a few remaining narrow gaps in 
absorption) in the wavelength interval ~ 7550 — 8100 A, 
between the patch of the Lya forest in the wavelength 
range ~ 6900 — 7550 A and the quasar signal at > 8100 
A, is quite str iking. (The Mg II doublet at A « 9180 A 



was noted by Fan et al. 2000; using a weighted average 



of several lines, we measure the absorber redshift to be 
z = 2.27865.) 

The sharp Ly/3 drop at A ~ 6910 A allows us to esti- 



1240 line as Lya. This unfortunately precludes the study 
of the red wing of the Lya absorption, which may contain 
useful information about the structure of the reionization 
front flMadau fc Rees 2000| , |Loeb fc Barkana 200 1| , [Barkana 
& Loeb 2001). Our redshift is supported by the possible 
Ly7 line and a drop near A w 6540 A, and a Lyman limit 
(clearly seen in the 2-d spectra images) at A « 6135 A. 
The effects of the blue BAL wing of [N V] and Lya may 
extend as far as A ~ 7900 — 8000 A, but this is not critical 
for our discussion below. 

3. DISCUSSION AND CONCLUSIONS 

There has been much recent progress in theoretical un- 
derstanding and mod eling of the re i onization era (excellent 
reviews include, e.g., Madau 200C, Loeb fc Barkana 2001 



Barkana & Loeb 2001, Shapiro 2001, etc.). A good un- 



derstanding of the structure and extent of the reionization 
is important by itself, as it reflects the earliest phases of 
structure formation, and also for the modeling of CMBR 
foregrounds at high angular frequencies. 

A simple picture of a clean-cut Gunn-Peterson trough 
now appears unlikely. The key issue is the dumpiness of 
the IGM, and the gradual development and clumpy dis- 
tribution of the first ion izing sources, either protoga laxies 
or early AGN (see, e.g., |Miralda-Escude et al. 2000| ). The 
reionization is expect ed to occur gradually as the UV cmis - 
sivity increases (cf. McDonald & Miralda-Escude 2001), 



and ionization overcomes the recombination rate, with the 
lowest density regions becoming fully reionized first. This 
is also sugges t ed by modern numerical simulations (e.g. 



Gncdin 2000, Ciardi et al. 2000, Umemura et al. 2001 



etcT) which predict an extended period of reionization, 
ranging from z ~ 15 to z ~ 5 or so. 

As we approach the reionization era from the lower 
rcdshifts, the Lya forest thickens, with an occasional 
transmission gap due to the intersection of ionized bub- 
bles along the line of sight; eventually a complete Gunn- 
Peterson trough is reached. In other words, the inherent 
non-uniformity of galaxy and structure formation is re- 
flected in the structure of the IGM phase transition corre- 
sponding to the reionization. Further complications arise 
from the proximity effect due to the source used to probe 
the IGM, i.e., a luminous quasar, and the nature, luminos- 
ity, and duration of other sources near the line of sight. 

This general picture is illustrated well in Fig. 20 of 
Loeb fc Barkana 2001 (which is t he same a s Fig. 40 of 
Barkana fc Loeb 2001| , or Fig. 6 of |Loeb 1999] ). The quali- 
tative correspondence with the observed spectrum of SDSS 
1044-0125 (Fig. 2) is striking. The dramatic increase in 
the opacity of the Lya forest at A > 7550 A, i.e., z > 5.2, 
is exactly what is expected in the approach to (or the tail 
end of) the reionization era. This is further illustrated 
in Fig. 3, which shows a dramatic thickening of the Lya 
forest absorption at these redshifts. A slightly different 
interpretation is that we are seeing a somewhat later evo- 
lutionary stage of the reionization process, i.e., remaining 
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islands of (mostly) neutral gas embedded in a growing sea 
of ionized hydrogen. 

The strong observed Ly/3 break may also be due in part 
to the patches of diffuse absorption we see so clearly in 
Lya. The overall appearanc e of the spectrum is sugges- 
tive of some of the models by Haiman & Loeb 1999, for the 
reionization redshift (in their terminology) a few percent 
lower than the source redshift. 

The dark portions of our spectrum at z ~ 5.2 — 5.6 have 
the flux consistent with zero, to within the photon noise. 
The lower limit to the optical depth (r.m.s., per pixel) is 
r > 4; if we average the flux over the dark portions of the 
spectrum, this limit is considerably higher, r > 6 or 7, de- 
pending on the redshift window used. Even if we assume a 
very conservative systematic sky subtraction error of ~ 1% 
of the continuum, the implied optical depth limit would be 
t > 4.6. The extrapolation of t he empirical sca ling laws 
found by Press et al. 1993 and Kim et al. 2001 to these 
redshifts suggests r ~ 2. This again indicates that we 
are seeing more absorption than would be expected from 
a simple extrapolation of the Lya forest. 

The few remaining transmission spikes are naturally in- 
terpreted as being due to the as-yet unpercolated reion- 
ization bubbles along the line of sight. An issue arises of 
whether the damping wings of the remaining neutral hy- 
drogen clouds would s uppress such transmission spikes (cf. 
Miralda-Escudc 1998). However, a s shown by Madau & 
Rees 100*0 and Cen & Haiman 200C , this depends strongly 
on the extent of the Stromgrcn spheres produced by the 
ionized sources, i.e., their luminosities and lifetimes, akin 
to the usual quasar proximity effect. Indeed, one expects 
some clustering of the first luminous sources, which are 
expected to form at the highest pe aks of the density field, 
due to biasing (cf. Djorgovski 1999] and references therein). 
The probably massive host of SDSS 1044-0125 is likely to 
have some luminous neighbors. 

We also checked whether any of the dark regions we see 
at z ~ 5.2 — 5.6 may be DLA systems with associated 
metallic lines, mainly the C IV doublet. None were found, 
as is expected from as yet unenriched gas. 

In their discovery paper, Fan et al. 2000 addressed the 
issue of reionization, and concluded that it is not detected 
in their data, on the basis of the few remaining transmis- 
sion spikes. Their were heavily binned spectrum seems 
roughly comparable to our LRIS data shown in Fig. 1, 
from which we cannot conclude much; a higher resolution 



spectrum, such as our spectrum shown in Fig. 2, is neces- 
sary. Furthermore, in their original redshift interpretation 
(which was due to the t hen unknown BAL nature of the 
object) Fan et al. 2000 may have mistaken the leftover 
quasar flux around the Lya as being a part of the Lya 
forest. 



After this paper was submitted, Becker et al. 2001 pre- 
sented spectroscopy of this and t hree additional q uasars 
at z ~ 5.82 - 6.28, discovered by |Fan et al. 200l[ They 
present a compelling evidence for a Gunn-Petcrson trough 
in the spectrum of the most distant quasar, suggesting the 
reionization epoch at z ~ 6, as anticipated here. This is a 
crucial result. 

However, the d iscussion of the remaining spectra by 
Becker et al. 2001 was limited by the available data, with 
relatively short exposure times. In order to increase the 
apparent S/N ratio, they binned the spectra by a factor 
of ~ 25 in wavelength, from ~ 0.154 A pixel -1 to 4 A 
pixel" 1 . The resulting loss of resolution makes it hard 
to detect dark windows similar to those seen in our ESI 
spectrum, especially given the differences in the S/N ra- 
tio. They further evaluated the mean optical depth in 
very wide bins, with Az — 0.2, which clearly precludes 
the detection of any dark windows with a smaller extent 
in redshift. Given these differences in the data and the 
analysis, we see no inconsistencies with our results. 

Taken together, the data so far suggest an extended and 
patchy end to the reionization era, as expected from mod- 
ern models of structure formation and of the reionization 
itself. Probing along more lines of sight with a high-S/N, 
high-resolution spectroscopy is necessary in order to place 
more quantitative observational constraints. Further in- 
sights will be obtained from direct detections of luminous 
sources responsible for the reionization at these and higher 
redshifts, their luminosity function and clustering proper- 
ties. 

We want to thank the staff of the W. M. Keck Observa- 
tory for their expert assistance. The LRIS data were ob- 
tained in the course of a collaborative project with F. Har- 
rison and P. Mao. We thank numerous colleagues whose 
constructive comments helped us improve the discussion 
presented in the paper. SGD acknowledges partial fund- 
ing from the Bressler Foundation. The work of DS was 
carried out at the Jet Propulsion Laboratory, Caltech, un- 
der a contract with NASA. 
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Fig. 1. — Spectrum of SDSS 1044-0125 obtained with LRIS (top), and the corresponding night sky (bottom). The three dashed lines 
represent a plausible range of the unabsorbed quasar power-law continua. 
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Fig. 2. — Spectrum of SDSS 1044-0125 obtained with ESI. Lyman series breaks corresponding to redshift z = 5.73 are indicated. In this 
rcdshift interpretation, most of the quasar's Lya line is absorbed. The dramatic change in the density of the Lyo forest at A J5 7550 A, and 
probably the corresponding Ly/3 absorption at A <J 6900 A, arc suggestive of the onset of the reionization era at z > 5.2. 
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Fig. 3. — A schematic illustration of the evolution of the Lyoj forest, in four redshift intervals. The top two windows are from an ESI 
spectrum of SDSS 1737+5828 at z = 4.94; the bottom two are from the spectrum of SDSS 1044-0125, presented here. The spectra have been 
renormaliscd by the best estimate of the continuum (for the bottom two panels, we used the middle power-law shown in Fig. 1). Almost all 
narrow spikes in the bottom panel are due to imperfect night sky emission lines subtraction. 



